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ABSTRACT 

We present our long term Chandra X-ray monitoring data for the gravitationally lensed quasar Q2237+0305 
with 20 epochs spanning 10 years. We easily detect microlensing variability between the images in the full 
(0.2-8 keV), soft (0.2-2 keV), and hard (2-8 keV) bands at very high confidence. We also detect, for the first 
time, chromatic microlensing differences between the soft and hard X-ray bands. The hard X-ray band is more 
strongly microlensed than the soft band, suggesting that the corona above the accretion disk thought to generate 
the X-rays has a non-uniform electron distribution, in which the hotter and more energetic electrons occupy 
more compact regions surrounding the black holes. Both the hard and soft X-ray bands are more strongly 
microlensed than the optical (restframe UV) emission, indicating that the X-ray emission is more compact than 
the optical, confirming the microlensing results from other lenses. 

Subject headings: accretion, accretion disks — black hole physics — gravitational lensing — quasars: individ- 
ual (Q2237+0305) 



1. INTRODUCTION 

Quasar microlensing is an important new tool for studies 
of the structure of quasar accretion disks based on the depen- 
dence of the microlensing magnification on the source size 
(e.g., Wambsganss 2006; Kochanek et al. 2007). Microlens- 
ing was first detected by Irwin et al. (1989) in Q2237+0305 
(Huchra et al. 1985). Afterwards, microlensing effects have 
been frequently reported in this system in different energy 
bands (e.g., Wozniak et al. 2000; Dai et al. 2003; Goicoechea 
et al. 2003; Vakulik et al. 2004; Anguita et al. 2008; Eigen- 
brod et al. 2008; Mosquera et al. 2009; O'Dowd et al. 2010). 
The quasar microlensing technique was first used to crudely 
determine emission sizes relative to the Einstein ring size (e.g, 
Wambsganss et al. 1990; Rauch & Blandford 1991; Falco et 
al. 1996; Mediavilla et al. 1998; Agol et al. 2000). However, 
since the magnification diverges on the caustics, the technique 
can be used to resolve emission regions arbitrarily smaller 
than the Einstein ring size, which allowed crude measure- 
ments of the X-ray continuum and FeKa emission regions 
(e.g., Chartas et al. 2002, 2007; Dai et al. 2003; Ota et al. 
2006; Blackburne et al. 2006; Pooley et al. 2006). Quasar 
microlensing models has been developed by different groups 
to interpret the microlensing results, where a major break- 
through occurred in Kochanek (2004) using a Bayesian Monte 
Carlo model. This approach and its derivatives have suc- 
cessfully reproduced microlensing light curves in the optical 
and X-ray bands, determined the emission sizes quantitatively 
(Morgan et al. 2008; Chartas et al. 2009; Dai et al. 2010; Me- 
diaville et al. 201 1), the stellar/dark matter fractions (e.g., Me- 
diavilla et al. 2009; Poindexter & Kochanek 2010a; Bate et 
al. 2011), and even the inclination angles of accretion disks 
(Poindexter & Kochanek 20 1 0b) . 

The X-ray continuum emission of quasars is generally be- 
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lieved to be created by the un-saturated inverse Compton scat- 
tering emission between soft UV photons from the disk and 
hot electrons in the corona, leading to the observed power- 
law X-ray spectra (e.g., Reynolds & Nowak 2003). Thus, the 
size of the X-ray emission traces the extent of the corona. Be- 
cause of its shorter variability time scales, the X-ray contin- 
uum emission region is believed to be smaller than the optical 
emission region. However, the simple variability argument is 
not conclusive as the variability could be dominated by lo- 
calized flares rather than disk-scale variability. Recent mi- 
crolensing results, however, have conclusively shown that the 
X-ray emission regions are an order of magnitude smaller than 
the (observed frame) optical emission regions (e.g., Morgan et 
al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et 
al. 2011), excluding models where the X-ray coronae covers 
large portions of the accretion disks. 

Besides the continuum emission, there are additional X-ray 
components such as the reflection hump, metal emission lines, 
soft X-ray excess, and jet emission, which may also contribute 
to the X-ray spectra of quasars. Quasar microlensing is one 
of the few tools that can determine the physical size of these 
different components, and microlensing effects have been de- 
tected in the FeKa emission line (e.g., Chartas et al. 2002, 
2007; Dai et al. 2003; Ota et al. 2006). In addition, if the 
hot electron distribution is not uniform, so we should observe 
chromatic, energy dependent microlensing of the X-ray con- 
tinuum. In this letter, we report results of our long term Chan- 
dra X-ray monitoring of Q2237+0305, and the first detection 
of chromatic microlensing of the X-ray continuum. 

2. OBSERVATIONS AND DATA REDUCTION 

In various combinations, we have now observed 
Q2237+0305 with ACIS (Garmire et al. 2002) on board 
the Chandra X-ray Observatory (Weisskopf et al. 2002) for 
20 epochs between September 2000 and November 2010 
(Table [l}. We reprocessed all these data using the CIAO 
4.3 software tools by removing the pixel randomization and 
applying a sub-pixel algorithm for event positions (Tsunemi 
et al. 2001). We filtered all events using the standard ASCA 
grades of 0, 2, 3, 4 and 6 with energies between 0.2 and 
8 keV (observer's frame), and further separated the events 
into soft (0.2-2.0 keV) and hard (2.0-8.0 keV) bands for 
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FIG . 1 . — Combined X-ray image of Q2237+0305 from 20 epochs observed 
between 2000 and 2010. The total integration time is 292 ks. 

our subsequent analyses. For each band, we extracted the 
total counts from all 4 images within 3" radius circular 
regions and the background counts from concentric annular 
regions with inner/outer radii of 10" and 20", respectively. 
We then used the PSF fitting method to model the relative 
count rates between the images, fixing the relative image 
positions to the HST positiono We binned the X-ray images 
in 0".0984, and fit these binned images using Sherpa to 
minimize the Cash C-statistic between the observed and 
model images. The final count rates were normalized by 
partitioning the background-subtracted total source counts 
using the (unnormalized) relative count rates for each image 
obtained from the PSF fits. We also obtained a stacked 
image of Q2237+0305 (Fig. 1) based on the image positions 
obtained from the PSF fitting results. 

We extracted spectra of the individual images in 0".5 ra- 
dius circles centered on the PSF fitting positions for each 
observation. We obtained combined spectra for each im- 
age by stacking the spectra and the corresponding rmf and 
arf files weighted by their exposure times, and analyzed 
the combined spectra using XSPEC V12. We fit the spectra 
of the individual (stacked) images using a simple power-law 
modified by Galactic absorption with Nh = 5.5 x 10 20 cm~ 2 
(Dickey & Lockman 1990) and neutral absorption at the lens 
(zi = 0.0395) for the 0.4-8.0 keV band. We detect differen- 
tial absorption between the images at the lens redshift. After 
obtaining the best fits, we calculated the ratios of the photon 
count rates between the absorbed and unabsorbed models for 
the full, soft, and hard bands, respectively. These mean ra- 
tios were then applied to the count rates of individual epochs 
and images to estimate the absorption corrected photon count 
rates. The absorption corrected count rates are presented in 
Table [T] Zimmer et al. (2011) computed the photon count 
rates of Q2237+0305 in the 0.5-10 keV band for 10 Chandra 
observations between 2006 and 2007; however, the authors 
did not correct for the absorption effects or divide them into 
energy bands. To double check our analyses, we also mea- 
sured the absorbed count rates for these observations in the 
0.5-10 keV band, and our results are consistent. 

3. ANALYSIS 
5 http://cfa-www.harvard.edu/glensdata/ 



TABLE 1 

Chandra Count Rates of Q2237+0305 Corrected for 
Absorption 
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NOTE. — Table Lis published in its entirety in the electronic edition, and a portion 
is shown here for guidance regarding its form and content. The full, soft, and hard bands 
are the observed energy bands between 0.2-8, 0.2-2, and 2-8 keV, respectively. 



Fig. [2] shows the time evolution of the flux ratios between 
various images. Since the time delays between different im- 
ages of Q2237+0305 are short, less than 1 day, from both 
theoretical and observational constraints (e.g., Schmidt et al. 
1998; Dai et al. 2003), the changes in the flux ratios are es- 
sentially unaffected by intrinsic variability, and thus must be 
due to microlensing by stars in the foreground lensing galaxy. 

We first fit the full, soft, and hard band light curves in Fig|2] 
to constants to test the significance of the microlensing signal. 
For the full band, the smallest \ 2 for no microlensing is for the 
case B/A, with a x 2 = 97.2 for 19 degrees of freedom (dof) 
and a null hypothesis probability of NP = 2 x 10~ 12 , while 
the largest is for the case of C/D with x 2 = 396.8 and NP = 
2 x 10" 72 . For the soft band, the smallest \ 2 = 70.2 is for B/A 
with NP = 8 x 10" 8 and the largest x 2 = 3 12.0 is for C/A with 
NP = 7 x 10" 55 . For the hard band, the smallest is x 2 = 44.35 
for B/D with NP = 0.00084, and the largest X 2 = 147.95 is 
for D/A with NP = 5 x 10" 22 . The null hypothesis is ruled out 
for all 6 image pairs for the full, soft, and hard bands. X-ray 
microlensing is detected with very high confidence. 

Next, we compare the hard and soft bands to test for chro- 
matic microlensing. Assuming no chromatic microlensing, 
we compute the x 2 f° r the hard and soft bands showing the 
same flux ratio evolution for three ways of combining the im- 
ages into pairs, i.e., B/A plus C/D, C/A plus B/D, and D/A 
plus C/B. The x 2 value for B/A+C/D is X 2 = 89.8 for 40 dof 
with a null hypothesis probability of NP = 1.1 x 10~ 5 . For 
C/A+B/D, we find x 2 = 91.0 with NP = 7.6 x 10" 6 . For the 
D/A+C/B pairs, we find x 2 = 50.5 with NP = 0.12. These 
results indicate that there is significant chromatic X-ray mi- 
crolensing in Q2237+0305. Fig. [2] also shows that the hard 
band microlensing variability of Q2237+0305 has a larger 
amplitude than that in the soft band. To compare the ampli- 
tudes of variability in the soft and hard bands, we calculate the 
r.m.s. variability for each image pair and show the results in 
Fig [3] The soft and hard microlensing variability amplitudes 
are different at the 99% confidence for the B/A plus C/D pairs, 
and at the 99.7% confidence for the C/A plus B/D pairs. In the 
last cases, D/A and C/B, the difference in the variability am- 
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FIG. 2. — Chromatic X-ray microlensing light curves of Q2237+0305 in the full X-ray band (2.0-8.0 keV, red circles), soft X-ray band (0.2-2 keV, green 
triangles), hard X-ray band (2—8 keV, blue diamonds), and OGLE V band (Wozniak et al. 2000; Udalski et al. 2006) and our SMARTS data (orange circles). 
Insets in the third column are blow-ups of the last X-ray data points (three observations made within 4 days). We slightly offset the observation dates for the soft 
and hard bands to make the difference more visible. 
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FIG . 3 . — X-ray microlening variability amplitudes in the full, soft and hard 
bands. 

plitudes are not well resolved. 

4. DISCUSSION 

As in previous studies, we detect significant microlensing 
in Q2237+0305 in the full X-ray band, where the flux ratio 
changes are dramatically different from those in the optical 
microlensing light curves, as shown in Fig. 2. These results 
confirm the earlier studies in other gravitational lenses that 
the X-ray continuum emission regions are more compact than 
the optical emission regions (Morgan et al. 2008; Chartas et 
al. 2009; Dai et al. 2010). By dividing the X-ray energy band 
(0.2-8.0 keV) into soft and hard bands, we detect for the first 
time chromatic microlensing of the X-ray continuum in the 
sense the hard band is more strongly microlensed than the 
soft. For our radio-quiet quasar at z s = 1 -69, the observed 0.2- 
8 keV band probes the rest frame 0.5-22 keV band. Our spec- 
tral analysis does not reveal a reflection or a soft X-ray excess 
component. The metal emission lines, e.g., Fe Ka, will fall in 



this band; however, their equivalent widths are usually small 
enough such that the broad band flux is little affected. There- 
fore, the chromatic X-ray microlensing occurs in the X-ray 
continuum emission in Q2237+0305. We note that if addi- 
tional X-ray emission components exist in the quasar spectra, 
the interpretation of chromatic X-ray microlensing could be 
more complicated. 

If the X-ray continuum emission was generated through the 
inverse Compton scattering of the soft UV photons from the 
accretion disk by hot electrons from the corona (Reynolds & 
Nowak 2003), our results indicate that the distribution of hot 
electrons in the corona is not uniform, and that the hot elec- 
trons are more compact than the cool electrons. Thus, even 
if the observed X-ray continua emission can be modeled as 
simple power-laws, our results indicate that the real emission 
mechanism has spatial structures. The evidence of structure 
in the "corona" favors complicated geometrical configuration 
such as the light bending (Fabian et al. 2005) or aborted jet 
models (Ghisellini et al. 2004). Previous microlensing analy- 
ses have constrained the size of the X-ray continuum to about 
~ I0r g (e.g., Morgan et al. 2008; Chartas et al. 2009; Dai et 
al. 2010), and our results in this paper suggest even smaller 
hard emission regions tracking the event horizon of the black 
holes. A detailed microlensing model of the microlensed X- 
ray emission from Q2237+0305 will be presented in a forth- 
coming paper. 
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